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We are computing the modifications for the scalar and pseudoscalar meson masses and mixing 
angles due to the proper accounting of fermionic vacuum fluctuation in the framework of generalized 
2 + 1 flavor quark-meson model and Polyakov loop augmented Quark Meson model (PQM). The 
renormalized contribution of the divergent fermionic vacuum fluctuation at one loop level, makes 
these models effective quantum chromodynamics (QCD)-like models. It has been explicitly shown 
that analytical expressions for the model parameters, meson masses and mixing angles, do not 
depend on any arbitrary renormalization scale. We have investigated how the incorporation of 
fermionic vacuum fluctuation in quark meson and PQM models qualitatively and quantitatively 
affects the convergence in the masses of the chiral partners in pseudoscalar (tt, jy, 77', K) and scalar 
(cr, oo, /o, k) meson nonets as the temperature is varied on the reduced temperature scale. Chiral 
symmetry restoration trends emerging from the temperature variations of the meson masses of 
the chiral partners and mixing angles, have been identified and compared with in different model 
scenarios. Comparison of present results in the quark meson model with vacuum term (QMVT) and 
PQM model with vacuum term (PQMVT) with the already existing calculations in the bare 2 + 1 
quark meson and PQM models, show that the restoration of chiral symmetry becomes smoother due 
to the influence of the fermionic vacuum term. We find that the melting of the strange condensate 
registers a significant increase in the presence of fermionic vacuum term and its highest melting 
is found in PQMVT model. The role of Ua{^) anomaly in determining the isoscalar masses and 
mixing angles for the pseudoscalar (77 and 77') and scalar (a and /o) meson complex, has also 
been significantly modified due to the presence of fermionic vacuum term. The interplay of chiral 
symmetry restoration effects and the setting up of C/a(1) restoration trends has been shown to be 
significantly modified by the incorporation of fermionic vacuum fiuctuation in the effective potential 
of the quark meson and PQM models. 

PACS numbers: 12. 38. Aw, ll.30.Rd, 12.39.Fe, ll.lO.Wx 



I. INTRODUCTION 

The strong interaction theory predicts that normal 
hadronic matter goes through a phase transition and 
produces a coUective form of matter known as the Quark 
Gluon Plasma (QGP) under the extreme conditions of 
high temperature and/or density when the individual 
hadrons dissolve into their quark and gluon constituents 
[ll-[l]. Relativistic heavy ion collision experiments at 
RHIC (BNL), LHC (CERN) and the future CBM 
experiments at the FAIR facihty (GSTDarmstadt) aim to 
create and study such a collective state of matter. Study 
of the different aspects of this phase transition, is a tough 
and challenging task because Quantum Chromodynamics 
(QCD) which is the theory of strong interaction, becomes 
nonperturbative in the low energy limit. However 
the QCD vacuum reveals itself through the process of 
spontaneous chiral symmetry breaking and phenomenon 
of color confinement. 

In the zero quark mass limit, chiral condensate works 
as an order parameter for the spontaneous breakdown of 
the chiral symmetry in the low energy hadronic vacuum 
of the QCD. For the infinitely heavy quarks, in the 
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pure gauge SUc{i) QCD, the Z(3) (Center symmetry 
of the QCD color gauge group) symmetry, which is 
the symmetry of hadronic vacuum, gets spontaneously 
broken in the high temperature/density regime of 
QGP. Here the expectation value of the Wilson line 
(Polyakov loop) is related to the free energy of a static 
color charge, hence it serves as the order parameter 
of the confinement-deconfinement phase transition Q. 
Even though the center symmetry is always broken 
with the inclusion of dynamical quarks in the system, 
one can regard the Polyakov loop as an approximate 
order parameter because it is a good indicator of the 
confinement-deconfinement transition 0] ■ 

The lattice QCD calculations (see e.g. [sj-flGj) 
give us important information and insights 
regarding various aspects of the transition, like the 
restoration of chiral symmetry in QCD, order of the 
confinement-deconfinement phase transition, richness 
of the QCD phase structure and mapping of the phase 
diagram. Since lattice calculations are technically 
involved and various issues are not conclusively settled 
within the lattice community, one resorts to the 
calculations within the ambit of phenomenological 
models [l7l - l26| developed in terms of effective degrees 
of freedom. These model investigations complement the 
lattice simulation studies and give much needed insight 
about the regions of phase diagram inaccessible to lattice 
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simulations. Lot of current effective model building 
activity, is centered around combining the features 
of spontaneous breakdown of both chiral symmetry 
as well as the center Z{3) symmetry of QCD in one 
single model (see for example j27l - l49| ). In these models 
chiral condensate and Polyakov loop are simultaneously 
coupled to the quark degrees of freedom. 

The behavior patterns of mesons and their properties 
in the hot and dense medium, have been investigated 
in the several two and three flavor NJL, PNJL models 
(e.g. js^-Ull) and also in the SU{2) version of linear 
sigma model (e.g. [HsI - fSTj ). Since the parity doubhng 
of mesons signals the restoration of chiral symmetry, 
these studies look for the emergence of mass convergence 
patterns in the masses of the chiral partners in 
pseudoscalar (tt, 77, 77', K) and scalar mesons (cr, ao, 
/o, k). We know that the basic QCD lagrangian has 
the global SUb.+l{3) x SUR-L{i) x Ua{1) symmetry. 
For the SU{3) Linear Sigma Model, several explicit 
as well as spontaneous symmetry breaking patterns of 
SUv{S) X S'C/yi(3), have been discussed by Lenaghan et 
al.in Rcf . j58| . Enlargin g th e Linear Sigma Model with 
the inclusion of quarks (59[ in the 2+1 flavor breaking 
scenario, Schaefer et al. studied the consequences 
of SU (3) chiral symmetry restoration for scalar and 
pseudoscalar meson masses and mixing angles, in the 
presence as well as the absence of C/yi(l) axial symmetry, 
as the temperature is increased through the phase 
transition temperature. It was shown by 't Hooft [6^ 
that the ?7a(1) axial symmetry does not exist at the 
quantum level and the instanton effects explicitly break it 
to ZA{Nf). Due to the Ua{^) anomaly, the 77' meson does 
not remain massless Goldstone boson in the chiral limit of 
zero quark masses and it acquires a mass of about 1 GeV. 
This happens due to the flavor mixing, a phenomenon 
that lifts the degeneracy between the tt and rj' which 
otherwise would have been degenerate with tt in [/ (3) 
even if the explicit chiral symmetry breaking is present. 
There is large violation in Okubo-Zwcig-Iizuka (OZI) 
rule for both pseudoscalar and scalar mesons and ideal 
mixing is not achieved because of strong flavor mixing 
between nonstrange and strange flavor components of 
the mesons [s^. Hence [^4(1) restoration will have 
important observable effects on scalar and pseudoscalar 
meson masses as well as the mixing angles. 

The effect of Polyakov loop potential on the behavior 
of meson masses and mixing angles has been studied by 
Costa et. al in the PNJL model j53| a nd by Contrera 
et.al in the nonlocal PNJL model 53] • Here in the 
NJL model based studies, mesons are generated by some 
prescription [s^ and the rj' is not a well defined quantity 
[6l|. It becomes unbound soon after the temperature is 
raised from zero. In the 2+1 flavor quark meson linear 
sigma model investigations by Schaefer et al. [H, [6^ . 
the mesons are the explicit degrees of freedom included 
in the lagrangian from the very outset and the ^^4(1) 
breaking 't Hooft coupling term is constant. Recently, we 
investigated the influence of the Polyakov loop potential 



on the meson mass and mixing angle calculations in 
the scalar and pseudoscalar sector, in the framework of 
generalized 2+1 flavor quark meson model enlarged with 
the inclusion of Polyakov loop [2^ [63l - l65j . 

The chiral symmetry breaking mechanism in the 
QM/PQM model is different from that of the NJL/PNJL 
model. In the NJL/PNJL model, the fermionic vacuum 
fluctuation leads to the dynamical breaking of the chiral 
symmetry. While in most of the QM/PQM model 
calculations, fermionic vacuum loop contribution to the 
grand potential got frequently neglected till recently [l^, 
Isil . [55ll56l. [59} because here, the spontaneous breaking of 
chiral symmetry is generated by the mesonic potential 
itself. Recently, Skokov et. al. incorporated the 
appropriately renormalized fermionic vacuum fluctuation 
j6a | in the thermodynamic potential of the two flavor 
QM model which becomes an effective QCD-like model 
because now it can reproduce the second order chiral 
phase transition at = as expected from the 
universality arguments [l7j for the two massless flavors of 
QCD. The fermionic vacuum correction and its influence 
has also been investigated in earlier works [67l - l70| . In a 
recent work [7l|, we generalized the proper accounting 
of renormalized fermionic vacuum fluctuation in the two 
flavor PQM model to the non-zero chemical potentials 
and found that the position of critical end point shifts 
to a significantly higher chemical potential in the fi and 
T plane of the phase diagram. Very recently, Schaefer 
et. al.jT^l estimated the size of critical region around the 
CEP in a three flavor PQM model in the presence of the 
fermionic vacuum term. Sandeep et al. also investigated 
the phase structure and made comparisons with lattice 
data in another recent 2+1 quark flavor study with the 
effect of fermionic vacuum term (73j . In a very recent 
work [t^ I , the present author explored and compared the 
details of criticality in the two flavor QM, PQM models in 
the presence and absence of fermionic vacuum correction. 

In the present work, the author will explore how the 
proper accounting fermionic vacuum correction in the 
QM and PQM models, qualitatively and quantitatively 
affects the convergence of the masses of chiral partners, 
when the parity doubling takes place as the temperature 
is increased through Tc and the partial restoration 
of chiral symmetry is achieved. We will also be 
studying the effect of fermionic vacuum correction on 
the interplay of SUa{S) chiral symmetry and C/yi(l) 
symmetry restoration in the presence as well as absence 
of Polyakov loop potential in QM model. Since we are 
lacking in the experimental information on the behavior 
of mass and mixing angle observables in the medium, 
a comparative study of these quantities in different 
model scenarios and circumstances becomes all the more 
desirable. 

The arrangement of this paper is as follows. In 
Secin] we recapitulate the formulation of the model. 
The description of grand potential in the mean held 
approach has been presented in Sec. IIIII We have 
derived the modification of meson masses due to the 
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fcrmionic vacuum correction and presented the scale 
independent calculation of the effective potential with 
the renormalized model parameters in Scc lIVI where 
the formulae for meson masses and mixing angles have 
also been discussed. In Sec|Vl we will be discussing 
the numerical results and plots for understanding and 
analyzing the effect of fcrmionic vacuum correction on the 
chiral symmetry restoration. Summary and conclusion is 
presented in the last Sec lVIl 



II. MODEL FORMULATION 

We will be working in the generalized three flavor 
Quark Meson Chiral Linear Sigma Model which has 
been combined with the Polyakov loop potential [28l . 
|63| - [65| . In this model, quarks coming in three flavor 
are coupled to the SUv{i) x SUAi'i) symmetric mesonic 
fields together with spatially constant temporal gauge 
field represented by Polyakov loop potential. Polyakov 
loop field <I>(x) is defined as the thermal expectation value 
of color trace of Wilson loop in temporal direction 



Nr 



= — Tr.Lt 



(1) 



where L(x) is a matrix in the fundamental 
representation of the SUc{S) color gauge group. 



L{x) = Vexp 



drAolx, t) 



(2) 



Here V is path ordering, Aq is the temporal vector field 
and /3 = [|. 

The model Lagrangian is written in terms of 
quarks, mesons, couplings and Polyakov loop potential 

z^($,$*,r). 



(3) 



where the Lagrangian in Quark Meson Chiral Sigma 
model 

^QM = <lf {iYDf_, - g Ta {(7 a + ijsT^a) ) Qf + C-m (4) 

The coupling of quarks with the uniform temporal 
background gauge field is effected by the following 
replacement = dp, ~ lA^, and A^, = (5^o^o (Polyakov 
gauge), where A^ = gsA1;^X°-/2. is the SUc{3) gauge 
coupling. Xa are Gell-Mann matrices in the color space, 
a runs from 1 • • • 8. qf — {u, d, s)^ denotes the quarks 
coming in three flavors and three colors, g is the flavor 
blind Yukawa coupling that couples the three flavor of 
quarks with nine mesons in the scalar (ctq, = 0^) and 
pseudoscalar {iTa, = 0~) sectors. 

The quarks have no intrinsic mass but become massive 
after spontaneous chiral symmetry breaking because of 
non vanishing vacuum expectation value of the chiral 



condensate. The mesonic part of the lagrangian has the 
following form 

Crn = Tr [dpM'^d^'M) - m^TT:{M'^ M) - Ai [Yy{M^ M)]^ 
-AaTr [AlHlf + c[det{M) + det{M^)\ 
+Tr [H{M + M^)] . (5) 

The chiral field M is a 3 x 3 complex matrix comprising of 
the nine scalars a a and the nine pseudoscalar tt^ mesons. 



M TaCa = Ta((Ta + ITT a) 



(6) 



Here Ta represent 9 generators of U{2>) with Tq = 

a = 0, 1 ... 8. \a are standard Gell-Mann matrices 



with Aq 



1. The generators follow U{i) algebra 



[Ta,Tb] = ifabcTc and {Ta.Tb} = dabcTc where fabc 
and dabc are standard antisymmetric and symmetric 
structure constants respectively with fabo = and dabo = 

1 (5ab and matrices are normalized as Tr (ToTb) = 

The S'C/l(3) X SUu{3) chiral symmetry is explicitly 
broken by the explicit symmetry breaking term 



H 



(7) 



Here H is a 3 x 3 matrix with nine external parameters. 
The ^ field picks up the nonzero vacuum expectation 
value, ^ due to the spontaneous breakdown of the chiral 
symmetry. Since ^ must have the quantum numbers of 
the vacuum, explicit breakdown of the chiral symmetry 
is only possible with three nonzero parameters /iq, and 
hs- We are neglecting isospin symmetry breaking hence 
we choose hg, /ig 7^ 0. This leads to the 2 + 1 flavor 
symmetry breaking scenario with nonzero condensates 
(To and (jg. 

Apart from ho and hg, the other parameters in the 
model are flve in number. These are the squared 
tree- level mass of the meson flelds m^, quartic coupling 
constants Ai and A2, a Yukawa coupling g and a 
cubic coupling constant c which models the C/a(1) axial 
anomaly of the QCD vacuum. 

Since it is broken by the quantum effects, the 
?7^(1) axial which otherwise is a symmetry of the 
classical lagrangian, becomes anomalous [zll and gives 
large mass to 77' meson (m^' = 940 MeV). In the 
absence of anomaly, 77' meson would have been 

the ninth pseudoscalar Goldstone boson, resulting due 
to the spontaneous break down of the chiral J7a(3) 
symmetry. The entire pseudoscalar nonet corresponding 
to spontaneously broken Ua{3), would consist of the 
three TT, four K, rj and ?/' mesons, which are the massless 
pure Goldstone modes when H = and they become 
pscudo Goldstone modes after acquiring finite mass due 
to nonzero H in different symmetry breaking scenarios. 
The particles coming from octet (oq, /o, k) and singlet 
(cr) representations of SUyi'i) group, constitute scalar 
nonet (cr, ag, /o, k). In order to study the chiral 
symmetry restoration at high temperatures, we will be 
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investigating the trend of convergence in the masses of 
chiral partners occurring in pseudoscalar (tt, 77, 77', K) 
and scalar (cr, ag, fo, nonets, in the 2 + 1 flavor 
symmetry breaking scenario. 



A. Polyakov loop potential 

The effective potential U ($, $*, T) is constructed such 
that it reproduces thermodynamics of pure glue theory 
on the lattice for temperatures upto about twice the 
deconfinement phase transition temperature. In this 
work, we are using logarithmic form of Polyakov loop 
effective potential [S^. The results produced by this 
potential are known to be fitted well to the lattice results. 
This potential is given by the following expression 



"'-(^''^^-n . _^<l,*<, + 5(r)ln[l-6$*<i> 



T4 



-4(<I>*3 + $3) - 3($*$)2] (8) 



where the temperature dependent coefficients are as 
follow 



where V is the three dimensional volume of the system, 
and /3 = y. For three quark flavors, in general, the 
three quark chemical potentials are different. In this 
work, we assume that S'C/y(2) symmetry is preserved and 
neglect the small difference in masses of u and d quarks. 
Thus the quark chemical potential for the u and d quarks 
become equal fix ~ fJ-u ~ l-id- The strange quark chemical 
potential is fiy ~ fig. Further we consider symmetric 
quark matter and net baryon number to be zero. 

Here, the partition function is evaluated in the 
mean-field approximation [2^ |5^, [H^, [63|. We replace 
meson fields by their expectation values (<&) Tq(7q + 
Tgag and neglect both thermal as well as quantum 
fluctuations of meson fields while quarks and antiquarks 
arc retained as quantum fields. Now following the 
standard procedure as given in Refs. [13, [H, |4^ [t^ 
one can obtain the expression of grand potential as 
sum of pure gauge field contribution U (<f>, $*, T). meson 
contribution and quark/antiquark contribution evaluated 
in the presence of Polyakov loop, 



f^MF(T,/i) = = C/((To,a8)+Z^($,$*,r) 



V 



(10) 



a{T) = ao + a, ('|] + ' b{T) = 63 ' ' 



The parameters of Eq. ([H]) are 



flo = 3.51 , 
02 = 15.2 , 



&3 



-2.47 
-1.75 



The critical temperature for deconfinement phase 
transition Tq = 270 MeV is fixed for pure gauge Yang 
Mills theory. In the presence of dynamical quarks Tq is 
directly linked to the mass-scale Aqcd, the parameter 
which has a flavor and chemical potential dep endence in 
fuU dynamical QCD and Tq ro(iV/, /i) 0,^. 



In order to study 2-1-1 flavor case, one performs following 
basis transformation of condensates and external fields 
from original singlet octet (0, 8) basis to nonstrange 
strange basis (x, y). 



V3 



'^0 - \/ 77^8 



(11) 

(12) 



Similar expressions exist for writing the external fields 
{hx, hy) in terms of {ho, hg). Thus the nonstrange and 
strange quark/antiquark decouple and the quark masses 
become 



III. GRAND POTENTIAL IN THE MEAN 
FIELD APPROACH 

We arc considering a spatially uniform system 
in thermal equilibrium at finite temperature T and 
quark chemical potential /x/(/ = u,d,s). The 
partition function is written as the path integral over 
quark/antiquark and meson fields [28l [soj 



Z = Trexp[-/3(H- ^ fifMj)] 



f—u,d.s 



J YlVaaViTaJ VqVqcxp ~ dr j d 



^QMS 



f—u,d.s 



(9) 



jTLx — Q — , m„ = q— 
y 2 ' y ^ ^ 



(13) 



Quarks become massive in symmetry broken phase 
because of non zero vacuum expectation values of the 
condensates. 

The mesonic potential in the nonstrange-strange basis 
reads. 



772 C 



2V2 



+ y'^>^^ + ^(2Ai + A2)a4 

+ i(2Ai+2A2)a^ , (14) 

The quark/antiquark Polyakov loop contribution is 
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Meson masses for Scalars and Pseudoscalars rj' ,rj,riNS,ris 




Pseudoscalar Meson masses 










+ Al(x2 + + (A2v-«-.)^2 v^c^ 


ml 


+ fa; + Jm^— 






Al(a;2 + y2) + (^2 _ ^ 2y2) 

— — 5mk-\- 


2 


^(7a" + 4y2a;j/ + 5?/") + (A2v - ^) (x' + y^) 




2 

Tip, 00 


Xi{x^ + 2/2) + + r ) + §(2a; + ^/2^/) 




- ^{V2x + y)- {x^Sm^ + y^Smy) 
^{5x^ - 4V2xy + 7y^) + (A2v - f) + 2y^) 
+ ^{y/2x - f ) - i^Sm^ + 2y^Smy) 






2 1/2c- , 2c- \ 

+m — -^(a; onix + y omy) 


2 

"Is, 88 




2 

"lp,88 


\i{x +y) + ^{x +4y ) - §(4x - V2i/) 








1 ™2 1 / 2 1 ^ 2 \ 




2|i(V2x2 -xy- V2y') + V2 (A^v - f ) (4 " v') 
+ 375 - V2y) - V2i^Sm^ - y^Srriy) 




"lp,08 


^(x-2-2y2)-§(x/2x-2y) 
_ ^^(^:2J^^ „ 2y'Smy) 




"^s(p),oocos^ Os^p) +m3(p)_gg sin^ O^^p) + 27Ti^(p)_o8 sin9s(p) cos 


?3(p) 


m2 


+ \i{x^ + y^) + - ^y 




"^s(p),ooSin^ ^.(p) +m^(p)_g8 cos^ 6l,(p) - 2m^(p)_og sin6',(p) cos 






+ Ai(x2 + y^) + - V2xy + 2y^) - §x 




H'2ml^pyoO + '^?(p),88 + 2^"i'(p),08) 




m™ 2 
"'p.OO 


+ Ai(x-2 + y2) + ^(a:2 + y2) + ^(2x + V^y) 




|('^s(p),oo + S'Tijfp) gg — 2\/2m^(p)_Qg) 




™m 2 

rripss 


+ Mx' + y^) + f (x^ + 4y2) - f ^ + 



TABLE I: The squared masses of scalar and pseudoscalar mesons appear in nonstrange strange basis. In this table x = a^, 
y = a„ A2V = A2. + A2+, n = f^, 5m. = 4nlog {jjj^), Sniy = 4nlog and 5m,± = ^^f;^ log (^) . 

The masses of nonstrange (Jns{vns), strange 0-3(773) mesons are given in the last two rows of second column. 



Model 


c[MeV] 


[MeV'^] 


Ai 


A2s 


[MeV^] 


hy [MeV''^] 


QM W/Ua{1) 


4807.84 


(342.52)2 


1.40 


46.48 


(120.73)^ 


(336.41)^ 


QMVT W/Ua{1) 


4807.84 


-(184.86)2 


-1.689 


46.48 


(120.73)^ 


(336.41)^ 


QM W/oUa{1) 





-(189.85)2 


-17.01 


82.47 


(120.73)^ 


(336.41)3 


QMVT W/oC/a(1) 





-(424.68)2 


-20.46 


82.47 


(120.73)^ 


(336.41)3 



TABLE IL parameters for rria — 600 MeV with and without (7a (1) axial anomaly term. 



written as, 



ovac 



-2iV 



(27r)3 
T[lng+ + In gj]} 



(15) 



The first term of the Eq. p5|) is the renormalized 
fermionic vacuum loop contribution. As evaluated in 
Ref . [H, [rij , it is given by 



gq 



\ M 



8^2 -/-vi^y- ^^^^ 

Here M denotes an arbitrary renormalization scale. In 
the second term and gJ are defined after taking trace 
over color space 



9f 



(17) 



(18) 



Ef ^ fj, and Ef is the flavor dependent single 



particle energy of quark/antiquark. 



(19) 



TOy is the flavor dependent quark mass. 



IV. MESON MASSES, RENORMALIZED 
MODEL PARAMETERS AND MIXING ANGLES 

In the effective potential Ea. ([TU|) . the pure mcsonic 
potential U((T^,crj,) and the renormalized fermionic 
vacuum term are relevant for fixing the vacuum (T = 







2 / 2 


2 2 / 
•niy^arriy^bl 


9" ml^ab/g'^ 




CTO 


1 „2 


2 
3 


1 „2 
3'^y 


1 

3 


CTl 


CTl 


12 
2^ 


1 








0-4 


0-4 





<y^ -\-\/2a'y 





\/2cr^4-2cry 




2aj;-<.^ 


erg 

O-Q 


crs 
0-8 


1 ^2 

Vlcr2 
6 ^ 


1 
3 

3 


2 ^2 
_v^cr2 

3 y 


2 
3 

V2 
3 


TTo 


TTo 





2 
3 





1 
3 


TTl 







1 








TV4 


7r4 





fTg. — \/2.l7y 











TTg 







1 

3 





2 
3 


vro 







3 
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TABLE III: First and second derivative of squared quark mass 
in nonstrange-strange basis with respect to meson fields are 
evaluated at minimum. Sum over two light flavors, denoted by 
symbol x, is in third and fourth column. The last two columns 
have only strange quark mass flavor denoted by symbol y. 
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0, /-i = 0) parameters m?, Ai, A2 and c. We write 



qq 



(20) 



The curvature of the grand potential in Eq. (|TO| at the 
global minimum gives scalar and pseudoscalar meson 
masses. 



(21) 



here subscript a = s, p; s stands for scalar and p stands 
for pseudoscalar meson and a, 6 = • • • 8. We note that 
the Polyakov loop decouples from the mesonic sector at 
T=0 and the curvature of the grand potential reduces 
to the second derivative of ft{(Tx,<Ty) in Eq. (|20l). In 
the PQMVT model, the meson mass modifications due 
to the fermionic vacuum contribution of the Dirac sea of 
quarks/antiquarks will be given by 



92 vac 




We write ^ 



and m^^^j, 



'// ^i,a,a as the first partial derivative 

dm^jr^^/d^a.b as the second partial derivative 
of the squared quark mass with respect to the meson 
fields ^a,b- These derivatives as originally evaluated 
in [5^ are given in Table IIIII The expressions for 
the vacuum meson masses and mixing angles which 
result purely from the second derivative of the mesonic 



potential \J{ax,crx, T, ^) are given in Table II of Ref. [28 1 
with their details discussed in Ref. (H,!!^. In the present 
calculation the vacuum masses that determine A2 and c 



are mt 



'"■pfiO 

that r. 



= m 
-Srn 

n2 



p,oo 
m 



and ?7i2_gg = 

i2 _ m 2 
r/' — '"■p,00 



K 



f (5?7ip_Qg^. We know 
^ The superscript 



m denotes the pure mesonic potential contributions to 
the pseudoscalar meson masses. We have collected the 
expressions of these masses in the last four rows of the 
third and fourth column of Table |lll After evaluating 



the mass modifications 5r 



''p.ab 



due to the influence of 



fermionic vacuum for a,b=0,l,4 and 8, we give a brief 
description of the evaluation of the parameters A2 and c 
in the following. 



rrif 



+ [^ + 2m} log 



/mf\ \ cl^mj 



(22) 



in 2 
777,7,- =777 



[m 2 I ni 2 




.T^X + 2V2y^Y) ; m^" = ml 



647r2 



x^X with X = ax,y = cTy 



1927r2 

A2 and c are obtained as 

3 (2/k - /.) m'°2 _ (2/^ + /,) ^^2 



{3x'X + 6y'Y) ; X = 1 + 41og (1^) , r = 1 + 41og (^-g^) (23) 



A2 



2(m-2^m-^)(/ 



K 



(3/2 + 8fK Uk - U)) [Ik - U) 



Jk - fn 



- A2 (2/k - ./.)(24) 



Putting the value of m™^, and 



A2s + A2+ + n 1 + 4 log 



^28 



9i2fK~U) 
2M 



in Eq.dMl) we find 



where A24 



!k [Ik - U) 



log 



2 Jk ~ fir 



327r2 



3 {2fK - U) ml ~ (2/k + h) ml - 2 (m^ + mg,) {Jk - U) 
(3f^ + 8fKifK-U))ifK-U) 



and c ~ 



mx — m 



Ik - fiT 



— -X2s{2fK-U) (25) 



r 



The fermionic vacuum modifications Sm"^ of the 
meson masses have a logarithmic dependence on the 
renormalization scale M. We denote the parameter 
A2 determined in the earlier QM and PQM model 
calculations in Ref. [M [H [111 by A2S. The proper 
renormalization of fermionic vacuum leads to an addition 
of a constant (A2++n) contribution to the A2S. Apart 



from it, the A2 in Eq. (|25p has a renormalization scale 
M dependent contribution in the present evaluation. 
Substituting this value of A2 in V{ax, cry), we find that the 
logarithmic M dependence of A2 completely cancels the 
scale dependence of term in Eg. ipU)) and the chiral 
part of the total effective potential becomes free of any 
renormalization scale dependence. We get the complete 
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cancellation of M dependence in the evaluation of c also 
and finally its value turns out to be the same as in the 
QM model. Further when this value of A2 is substituted 
in the expressions of meson masses ( as given in Table 
II of Ref. [2^) obtained from the second derivative of 
the pure mesonic potential U((Ta;, ctj^), the logarithmic 
M dependence of A2 neatly cancels with the scale 
dependence already existing in the mass modifications 
(5™^ ab^ '^^'^ fermionic vacuum correction and the 

final expression of meson masses becomes free of 

any scale dependence when these two contributions are 
added together. The expressions of the scale independent 
scalar and pseduscalar meson masses that we will be 
using in the present calculation are given in Table |lll 
+ A2+ has been denoted by Av. 

The chiral part of the PQM/PQMVT model has 
six input parameters and therefore require six known 
quantities as input. In general ttZtt, Wif, the pion and 
kaon decay constant Z^, /x, mass square of t]' and 
TOct are used to fix these parameters. The parameters 
are fitted such that in vacuum, the model produces 
observed pion mass 138 McV, kaon mass 496 McV and 
m'^ = 963(138) MeV, = 539(634.8) McV for the case 
with the presence (absence c=0) of axial anomaly term c. 
We have already obtained the A2S and c, numerical values 
of A2+ and n are obtained using /^r = 92.4, fx = 113 MeV 
and Nc = 3. Using the expression of m^, we express 
in terms of Ai. With = 600 MeV, the expression of 
TO^ given in Table HI] is exploited to yield the numerical 
value of Ai. In the present work, the A2S and c are 
the same as in [s^, the value of and hy are also 
not affected by the fermionic vacuum correction. The 
parameters which are modified by the fermionic vacuum 
correction are , Ai and A2. Table ITlAl summarizes the 
parameters in different model scenarios. 

We finally obtain the following scale independent chiral 
part of the effective potential as: 

f^K, CTy) = ^ (c^^ + <^l) - hxCTx - hy<7y - -^^^l^ V 
^1 ,2„2 , -^1 /'„4 , „4^ , (•^2v+") / 4 , r,^4\ 
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4 / 4 

no^ I erg 

' 2 °^[{2fK-U)J '^'^"'^^\{2fK-U) 



nat loR ( .^^"""^ , 1(26) 



Now the thermodynamic grand potential in the 
presence of appropriately renormalized fermionic vacuum 
contribution in the PQMVT model will be written as 



an 



on 

do-,, 



an 



an 







CT I, = (Tai , a = cr y . = $ , $ * = <I) * 

(28) 

The diagonalization of (0 - 8) component of mass 
matrix gives masses of a and /o mesons in scalar sector 
and masses of 77' and 77 in pseudoscalar sector. The scalar 
mixing angle 9s and pseudoscalar mixing angle 6p are 
given by, 



tan 26n 



2^0,08 



'0,88 



(29) 



The meson masses are further modified in medium at 
finite temperature by the quark contributions in the 
grand potential. In order to calculate the second 
derivative Eq.()2ip for evaluating the finite temperature 
quark contribution in the presence of Polyakov loop 
potential, the complete dependence of all scalar and 
pseudoscalar meson fields in Eq.® has to be taken into 
account. We have to diagonalize the resulting quark mass 
matrix. The expression for the meson mass modification 
due to the quark contribution at finite temperature in 
the QM and PQM model, has already been evaluated 
respectively in Ref. [s^ and [l^ In the following, we 
give the expression of mass modification due to the quark 
contribution at finite temperature in the PQM model as: 



3E 

f=x,v 



d^p 1 



2Ej 



-(B++B-) 



The notations and have the following 



definitions 



2EfT 



(30) 



^/ = = 

and B^ — 3{A^)^ — C^, where we again define 



(31) 



(32) 



nMF(r,Ai;a„(7y,$,$*) = U{T;'S>,<i>*)+n{a,,ay) + 

r!T(r,M;a„a„$,<i>*) (27) 



9f 



(33) 



One can get the quark condensates ax, cry and 
Polyakov loop expectation values $, $* by searching the 
global minima of the grand potential for a given value of 
temperature T and chemical potential fi. 



CZ 



q,*^-Ej IT ^ 4$g-2iJ^ IT ^ 3g-3iJ^. IT 



9f 



(34) 
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FIG. 1: The variations of nonstrange ax, strange Uy condensates with respect to the relative temperature scale (T/T^) at zero 
chemical potential {fi = 0) in the QM,QMVT ,PQM and PQMVT models have been shown. The magenta lines with dash 
double dots , dark green lines with thick long dash and red lines with dash dot, represent the respective variations in QM, 
QMVT and PQM models while the solid black lines represent the PQMVT model variations. Four such lines in the upper 
half of the figure represent the strange ay condensates while the same line types in the lower half represent the nonstrange ax 
condensates for computations done with the axial anomaly (non-zero c). In the upper half of the figure, the first lower line with 
dark blue solid circular dots and then the second lower line with dark green solid triangular dots show the respective variations 
of the ay in the PQMVT and PQM models while the upper line with dark red solid circular dots and the uppermost line with 
black solid triangular dots represent the respective ay variations in the QMVT and QM models when axial anomaly term is 
absent i.e. c — 0. The expectation value of the Polyakov loop {$), is shown in the right hand side plots, dark blue small dash 
line represents the $ variation in the PQMVT model while the double dash magenta line represents the $ variation in the 
PQM model. 



QM QMVT PQM PQMVT 

T^{MeV) 146.1 171.1 205.8 216.5 
Ti^(MeV) 248.3 ± .5 247.8 ± 1.0 274 ± 1.5 269. ± 1.5 
(MeV) - 205.6 205.6 



TABLE IV: The characteristic temperature (pseudo critical 
temperature) for the chiral transition in the nonstrange 
sector T^, strange sector T^- and confinement-deconfinement 
transition Tf , in QM, QMVT, PQM and PQMVT model. 



V. FERMIONIC VACUUM CORRECTION AND 
CHIRAL RESTORATION 

We are investigating the effect of fermionic vaeuum 
fluctuation on the restoration of chiral symmetry when 
it is properly accounted for in the 2 + 1 flavor quark 



meson model and PQM model at finite temperature 
and zero chemical potential with and without axial 
Ua{^) breaking. We have compared the results of 
present computations in the QMVT and PQMVT 
models with the already existing calculations in the 
quark meson model and PQM model [2^, [s^. The 
interplay of the effect of Ua{^) axial restoration 
and chiral symmetry restoration in the influence of 
fermionic vacuum fluctuation has been demonstrated 
and compared with in different model scenarios through 
the temperature variation of strange, nonstrange chiral 
condensates, meson masses and mixing angles. The 
C/^(l) axial breaking term is constant throughout the 
computation. The value of Yukawa coupling g has been 
fixed from the nonstrange constituent quark mass = 
300 MeV and is equal to 6.5. This predicts the strange 
quark mass i2 433 MeV. 
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A. Condensates and fermionic vacuum correction 

The solutions of the gap equations fEa l^5|) yield 
the temperature dependence of the Polyakov loop 
expectation value ($), nonstrange and strange 
condensates at zero chemical potential and the inflection 
points of these order parameters respectively give the 
characteristic temperature (pseudocritical - temperature) 
for the confinement - deconfinement transition T*, the 
chiral transition in the nonstrange and strange sector 
T^. Numerical values of pseudo critical temperatures 
for different transitions in quark meson (QM), QMVT, 
PQM and PQMVT models, are given in Table El It 
is evident from the Table |V] that the chiral crossover 
transition at /i = shifts to higher temperatures due 
to the fermionic vacuum correction. Since the thermal 
fermionic contributions to the effective potential are 
suppressed due to the presence of the Polyakov loop, the 
chiral transition gets delayed and as we already know, 
the ^ = crossover occurs at a significantly higher 
temperature (28l . [63j in PQM model. The incorporation 
of fermionic vacuum fluctuation in the effective potential 
of QM model also changes the model parameters such 
that the chiral transition gets delayed and the QMVT 
model increases by 25 MeV over its QM model 
value, only due to the effect of the fermionic vacuum 
correction. Further its influence, gives highest value 
of in PQMVT model. Here it is relevant to recall 
that larger values of nia and Tq also shift the chiral and 
deconfinement transition temperatures to higher values. 
We take m„ and Tq as 600.0 and 270.0 MeV respectively 
in the present calculation in order to compare results 
with our earlier work in the QM and PQM models in Ref. 
[28j . Further with these parameter sets, using different 
versions of Polyakov loop potential in PQM model, 
one gets coincident chiral and deconfinement crossover 
transitions [63j . We notice that the fermionic vacuum 
correction does not affect the confinement-deconfinement 
crossover transition temperature and T* has the same 
value of 205.6 MeV in both the models PQM as weh 
as PQMVT. But the chiral crossover transition in the 
non-strange sector does not remain coincident with the 
deconfinement transition due to the fermionic vacuum 
fiuctuation and registering an increase of 11 MeV over 
the deconfinement crossover temperature T*=205.6 
MeV, the PQMVT model becomes 216.5 MeV, as 
shown in Table [V] 

Wc compare the QMVT(PQMVT) model calculations 
with the corresponding results in the QM(PQM) 
model on a relative temperature scale T/T^. This 
is justified because absolute comparison of the 
characteristic temperatures between two models of the 
same universality class can not be made according 
to Ginsburg-Landau effective theory [53|. Since our 
investigation is focused on identifying the effective 
symmetry restoration trends due to the influence of 
fermionic vacuum fluctuations, we will be comparing the 
mesonic observables below and above T^. 



The condensates start with the values = 92.4 MeV 
and ay =94.5 MeV at T = in FigH] . The influence of 
Ua{^) anomaly (with zero or non-zero value of c) on the 
behavior of the nonstrange condensate is negligible. The 
ax variation which is sharpest in the T /T^ = 0.8 to 1.4 
range in the PQM model (dash dot line in red), becomes 
smoother in the PQMVT model (solid black line) due 
to the effect of fermionic vacuum term. In the absence 
of Polyakov loop, the QM model (dash double dot line 
in magenta) a^ temperature variation becomes quite 
smooth in the QMVT model (thick long dash line in dark 
green) , only because of the fermionic vacuum correction. 
The variation of strange condensate ay in PQM model 
(dash dot line in red) shows an early and signiflcant 
melting due to the influence of Polyakov loop potential 
when calculations are done with axial anomaly [28] . It is 
evident from the upper half black solid line in Fig[T] that 
the corresponding melting of ay is enhanced due to the 
fermionic vacuum correction and we obtain the largest 
degree of strange condensate melting in PQMVT model. 
In the influence of fermionic vacuum correction only, the 
modiflcation in the melting of strange condensate ay is 
quite larger when we compare the upper half thick long 
dash line (dark green) in QMVT model with the dash 
double dot line (magenta) obtained in QM model. When 
we compare the respective temperature variations of ay 
represented by the filled circular dots line in dark blue 
(red) color for the PQMVT (QMVT) model with the 
filled triangular dots line in dark green (black) color for 
the PQM (QM) model, we conclude that the melting of 
the strange condensate gets a little reduced in the same 
proportion in the absence of axial anomaly term i.e. for 
the c = case in aU the models PQMVT and PQM, 
QMVT and QM. 

Curves starting from the right end of the plot 
represent the temperature variation of the Polyakov loop 
expectation value ($). Line with small double dash 
(magenta) denote the $ variation in PQM model while 
line with small dash (dark blue) depict the PQMVT 
model $ variation. Here we recall that the improved 
ansatz of the logarithmic polyakov loop potential [1, 
[3^ . Issl . lisj avoids the $ expectation value higher than 
one and hence describes the dynamics of gluons more 
effectively. 

The chiral crossover transition in the strange sector 
is lot more smooth and weak in comparison to the 
corresponding chiral transition in non-strange sector in 
all the models due to the large constituent mass of the 
strange quark = 433 MeV. The variation of the 
temperature derivative of ay shows two peaks in all the 
models, the first peak is higher and sharper because it 
is driven by the chiral crossover transition dynamics in 
the non strange sector. The crossover temperature in 
the strange sector is identified by locating the position 
of the second peak which is quite broad, smooth and 
fiat over a small temperature range in all the models. 
Though we have not shown the temperature variation 
of the temperature derivative of ay in this work, we 
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FIG. 2: Mass variations for the chiral partners {a, n) and (ao, rj') on the reduced temperature {T/T^) scale at ^ = 0, are 
plotted in Fig 2(a) for the PQMVT and QMVT model and the corresponding mass variations in the PQM and QM model, are 
plotted in Fig 2(b) The axial anomaly coefficient c has a constant non-zero value in these computations. The a and ao mass 
variations are depicted by solid line magenta color plots and dark green line with dash dot plots respectively in the PQMVT 
(PQM) model and the QMVT (QM) model in the left panel (right panel) while the vr and rj' mass variations are denoted by 
the dash double dot black line plots and double dash red line plots respectively in the PQMVT (PQM) model and the QMVT 
(QM) model in the left panel (right panel). 
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FIG. 3: The line types in Fig |3(a) and in Fig 3(b) represent the same mass variations as depicted in Fig[2]but here in these 
computations the axial anomaly term is absent i.e. c = 0. 



show the ambiguity in the identification of second peak 
by indicating the flatness range for in Table |Vl As 
discussed above, the earliest and largest but smoother 
melting of strange condensate is obtained in the PQMVT 
model with T^= 269.0 ±1.5 MeV. It will have an 
interesting physical consequence in the early setting up 
of a smoother mass degeneration trend in masses of the 
chiral partners {K, k) and {rj, fo) and in the early 
emergence of a smoother ?7a(1) restoration trend. 



B. Meson Mass Variations 

The expressions of the vacuum masses in the QMVT 
model evaluated after including the contributions from 
the fermionic vacuum fluctuation for the scalar and 
pseudoscalar mesons, are given in Table |TT1 The 



corresponding vacuum meson mass expressions for the 
pure QM model (with no fermionic vacuum correction) 

When 
due to 
model 



28|. 



are already given in Table II of Ref. 
the finite temperature mass modifications 
the quark- ant iquark contribution in the QM 
(evaluated in [53|), arc added to the vacuum mass 
expressions in the QMVT and QM model, we get the 
total value of the finite temperature meson masses in the 
respective models. Similarly when the finite temperature 
mass modifications due to the qq contribution calculated 
in the presence of Polyakov loop potential in Eq. ([50]) , are 
added to the vacuum masses as obtained in the QMVT 
and QM model, we get the final value of the meson masses 
in the PQMVT and PQM model. In the presence of 
axial ?7a(1) anomaly, the respective mass variations in 
the PQMVT as well as QMVT model and the PQM 
as well as QM model for the chiral partners (cr,7r) and 
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{ao,T]'), are plotted in Fig 2(a) and Fig 2(b) respectively restoration trend 



and the corresponding mass variations for the chiral 
partners (?7,/o) and (K^n) are plotted respectively in 
Fig 4(a) and Fig 4(b) In the absence of the J7a(1) axial 
anomaly (i.e. for c=0), the corresponding model plots for 
the chiral partners ((T,7r) and (09,77') are given in Fig 3(a) 
and Fig |3(b)| respectively while Fig 5(a) and Fig 5(b) give 
the respective plots in the same model sets for the chiral 
partners (77, /o) and (K,k). In the influence of fermionic 
vacuum correction, the sharpest mass degeneration in 
PQM model computations in Fig 2(b) for {a,Tr) and 
(00,77') mesons, becomes quite smooth in the PQMVT 
model in Fig 2(a) and here the most smooth mass 
degeneration can be seen in the QMVT model due to 
the absence of Polyakov loop potential. We conclude 
from the behavior of these chiral partners that the 
net effect of the fermionic vacuum correction in the 
PQM and QM model, is to make a smoother but more 
effective occurrence of chiral SUl{2) x SUr{'2) symmetry 
restoration. Further in comparison to PQM and QM 
models, the degenerate a and tt meson masses in the 
PQMVT and QMVT models, show closer convergence 
towards the degenerate masses of og and 77' for higher 
T/TX > 1. 

We notice the similar trend of smoother but more 
effective mass degeneration for the chiral partners (7/,/o) 
and {K, n) in Fig 4(a)| due to the fermionic vacuum 
effect. In the PQM (QM) model, the /o mass intersects 
the r] mass nearly at T/T^ = 1.4(1.8) and becomes 
smaller than the developing a kink like structure after 
crossing it for higher val ues of the reduced temperature 
T/TX > 1.4(1.8) in Fig|4(b)l Further the /o meson 



becomes degenerate with (A", n) and -q again at T/T'^ > 
1.8(2.3). It is important to emphasize that this kink-like 
mass crossing behavior of /o meson altogether disappears 
from the PQMVT (QM VT) m odel mass variations of 
/o and 77 mesons in Fig 4(a) due to the noteworthy 
effect of the fermionic vacuum correction. Here the /g 
meson becomes degenerate with the {K, n) and r] mesons 
earlier at T/T^ > 1.6(1.9) and remains so forever. This 
trend of mass degeneration reflects the effect of fermionic 
vacuum fluctuation on the chiral symmetry restoration 
in the strange sector and it results due to the smoother 
but larger (largest in the PQMVT model) melting of 
the strange condensate in FigH] Here it is worthwhile 
to recall that the breaking of C/^(l) axial symmetry, 
generates the mass gap between the two sets of the 
chiral partners, (cr, tt) and (oq, 77') i.e. m-^ ~ m„ < 
"niaa — ''^T'ri' for T/Tq > 1. This results due to the 
opposite sign of the anomaly term {■s/2c(jy) in the scalar 
and pseudoscalar meson masses. Hence the mass gap 
reduction will be larger due to the larger melting of Uy 
for higher T/T^ > 1 and the mass of the degenerated cr, 
TT mesons will converge more closely to the higher mass of 
the already degenerated o q, 77' m esons. We have already 
seen this behavior in Fig 2(a) Thus the inclusion of 



Here we will be discussing the variations in the masses 
of the chiral partners when the Ua{^) axial anomaly 
is absent (c=0). Comparing the plots in Fig 3(a)| and 
Fig 5(a) respectively with those in Fig |3(b)| and Fig |5(b)[ 
we again find, the same smoother but more effective 
mass convergence trend for the mass degeneration that 
we identify as the effect generated by the inclusion of 
fermionic vacuum fluctuation in the PQMVT and QMVT 
model. In the absence of the anomaly, the = ttt,^ in 
the vacuum and it stays the same at all temperatures 



in Fig 3(a) and Fig 3(b) Further the mass gap between 
the chiral partners (cr, tt) and (ag, rj') becomes zero and 
all the four mesons become degenerate at T/T^ = 1.0 
in all the models . The T /T^ numerical value, where 
the if, K and 77 masses degenerate in different models, is 
not influenced by the Ua(X) anomaly as expected since 
the nonstrange condensate does not have any anomaly 
dependence. Further, in Fig |5(b)| the /o mass variation 
in the PQM (QM) model, does not show the crossing 
behavior similar to the corresponding mass variation 
computed in the presence of the axial anomaly. Here 
the /o mass does not become completely degenerate with 
the Tn,,, though it becomes very close (nearly touches) 
to the 77 mass variation when T jT^ ~ 1.6(2.0) and 
afterwards /o takes slightly larger value than m^. In 
the PQMVT (QMVT) model in Fig(5(a)l the /o mass 



fermionic vacuum fluctuation in the PQM (QM) model 
also effects an early and smoother set up of the Ua^X) 



variation becomes completely degenerate with the tti^ 
when T/T^ ~ 1.9(2.3) and does not show any crossing 
behavior. Further in the absence of anomaly we also 
note that the mass of /o in vacuum increases by about 
60 MeV. 

Here we mention another noteworthy flnding. In 
the influence of the fermionic vacuum correction, the 
scalar particle vacuum mass increases to 1086.26(917.93) 
MeV for the oq meson and decreases to 1143.92(1203.16) 
MeV for the /o meson in the presence (absence) of axial 
anomaly in the PQMVT and QMVT models from the 
respective vacuum mass value of rriao =1028.7(850.5) 
MeV and 77^/„ = 1221. 1(1282.3) MeV in the PQM and 
QM model. We recall further that the mass variations 
of scalar cr and /o show the kink around T/T^ = 1.8 in 
the QM model while it is seen around T/T'^ = 1.4 in 
PQM model. The kink generation in the PQM and QM 
model results because the meson masses interchange their 
identities for higher values on the reduced temperature 
scale [2^, • It is worth emphasizing that the kink in the 
scalar cr and /q meson masses, gets completely washed 
out by the presence of fermionic vacuum fluctuation in 
the PQMVT and QMVT model and we get a completely 
smooth mass degeneration trend in these models without 
any crossing or anti-crossing of meson masses. In order 
to have a proper perspective of the kink behavior in the 
PQM and QM model mass variations and of the complete 
washing out of the kinks in PQMVT and QMVT model 
results, one has to investigate, analyze and compare the 
scalar and pseudoscalar meson mixing angles. 
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FIG. 4: Mass variations for tlie chiral partners (77, /o) and (K, k) on the reduced temperature (T/T^) scale at /x = 0, are 
plotted in Fig 4(a) for the PQMVT and QMVT model and the corresponding mass variations in the PQM and QM model, 
are plotted in Fig |4(b)"[ The axial anomaly coefficient c has a constant non-zero value in these computations. The r; and fo 
mass variations are depicted by the solid line in magenta color plots and dash dot line dark green color plots respectively in 
the PQMVT (PQM) model and the QMVT (QM) model in the left panel (right panel) while the K and k mass variations are 
denoted by the dash double dot black line plots and double dash red line plots respectively in the PQMVT (PQM) model and 
the QMVT (QM) model in the left panel (right panel) 
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FIG. 5: The line types in Fig 5(a) and in Fig 5(b) represent the same mass variations as depicted in Fig|4]but here in these 



computations the axial anomaly term is absent i.e. c = 0. 



C. Meson Mixing Angle Variations 



We will finally be investigating the behavior of the 
scalar and pseudoscalar Qp mixing angles on the 
relative temperature scale. In the presence of axial 
anomaly, the lower(upper) solid lines in magenta color 
and dash dot black lines depict the 9p{0s) variations 
respectively in the PQMVT(P QM) and QMVT(QM) 

|6(b)[ ). In the 



ot line in dark 



model computations in the Fig 6(a) (Fig 
absence of axial anomaly, the dash double d( 
green color and double dash line in red color represent 
the scalar 9$ mixing angle variations respectively in the 
PQMVT( PQM ) and QMVT(QM) model computations 
in the Fig 6(a) (Fig |6(b)"] ). The pseudoscalar 9p mixing 
angle variations for this case are constant and are shown 
by black dash line and filled circular dots in blue color 



respectively in the PQMV T(PQ M) and QMVT(QM) 
model calculations in the Fig 6(a) (Fig |6(b)"| ). Comparing 
the PQMVT and QMVT model variations of 9p and 9s 



in Fig 6(a) with the corresponding PQM and QM model 
results shown in Fig |6(b)| we infer that the fermionic 
vacuum correction significantly modifies the axial Ua{'^) 
restoration pattern. 

The nonstrange and strange quark mixing is strong, 
and at T = one gets the pseudoscalar mixing 
angle 9p ~ —5° in all the models which remains 
almost constant in the chiral broken phase. The 
0P va riation near T/T^ = 1 in the PQMVT model 
in Fig 6(a) develops a small dip and then smoothly 
starts the approach toward the ideal mixing angle 
9p — > arctan^ ^ 35°, the corresponding $p = 
90°. Here <I>p is the pseudoscalar mixing angle in the 
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QMVT w U.(l) 
PQMVT w Ui(l) 
QMVT w/oUi(l) 
mMVT w/o Ujh) 
QMVT w/oUi(I) 
RjMVTw/oUjO) 



0.2 0.4 0.6 0.8 



1.2 1.4 1.6 



2 2.2 2.4 




,4 1.6 1.8 2 2.2 2.4 



(a)Scalar and Pscdoscalar mixing angle variations in PQMVT 
and QMVT model 



(b)Scalar and Psedoscalar mixing angle variations in PQM 
and QM model 



FIG. 6: In the presence of axial anomaly, the lower(upper) magenta color solid lines and dash dot black lines depict the 9p{6s) 
variations respectively in the PQMVT and QMVT model computations in t he Fi g f6(a)| while the same line types represent the 
corresponding variations respectively for the PQM and QM model in (Fig 6(b) I. In the absence of axial anomaly, the dash 
double dot line in dark green color and double dash line in red c olor represent the scalar 6s mixing angle variations respectively 
in the PQMVT and QMVT model computations in the Fig 6(a) while the same line types represent the corresponding variations 



respectively for the PQM and QM model in (Fig 6(b) I. The pseudoscalar Op mixing angle variations for this case are constant 



and are shown by black dash line and filled circular dots in blue color respectively in the PQMVT and QMVT model calculations 
in the Fig 6(a) while the same line types represent the corresponding variations respectively for the PQM and QM model in 
Fig|6(b)l 



strange nonstrange basis (see Ref.|59l| for details). In 
computations with the presence of axial anomaly for 
T/T^ > 1, the pseudoscalar mixing angle approaches its 
ideal value more smoothly in the PQMVT and QMVT 
model when compared with the corresponding result in 
the PQM and QM model in Fig |6(b)[ This approach is 
sharpest in the PQM model. 

The rj and r/ mesons become a purely strange rjs and 
nonstrange rji^s quark system as a consequence of the 
ideal pseudoscalar mixing which gets fully achieved at 
higher values of the reduced temperature. In order to 
show this and make comparisons, the mass variations 
for the physical 77, 77' and the nonstrangc-strangc ?77V5, 
77s complex, arc p lotted for the PQMVT and QMVT 
model in Fig |7(a) and t he pl ots for the PQM and QM 
model are shown in Fi g|7(b) Mass formula m^„e and 



'J/s 



are given in Table IIII In the to,j' approach to m 



Vns 



and the m,, approach to ruris around T/T^ = 1, the 
most smooth and smoother mass convergence trend is 
see n resp ectively in the QMVT and PQMVT model in 
Fig 7(a)| because of the influence of fermionic vacuum 
fluctuation. This mass convergence trend is sharp most 
in the PQM model in Fig |7(b)| due to the influence of 
Polyakov loop potential. The smallest mass difference 
between the mrj' (which is degenerate with m,-!^^) and 
rriri (which is degenerate with ) for T/T^ > 1 results 
in the PQMVT model. It means that we are getting the 
most effective J7a(1) restoration trend in the PQMVT 
model for the pseudoscalar sector. 

Due to the effect of fermionic vacuum correction, the 
scalar mixing angle 9$ value in vacuum (i.e. at T = 0) 
decreases to 11.98(14.75) degree in the presence(absence) 



of axial anomaly in the PQMVT and QMVT mod els in 
Fig 6(a) from its value of 19.86(21.5) degree in Fig |6(b)| 
for PQM and QM models. Growth of 9s to its ideal 
value starts near T/T^ = 1 in all the models but this 
growth is smoother in PQMVT and QMVT model. In 
the chirally symmetric phase of the PQMVT and QMVT 



model in Fig 6(a) the scalar mixing angle ^5 smoothly 
approaches the ideal value of {9s 35°) mixing angel for 
higher temperatures for both the cases with and without 
the axial anomaly. We thus conclude that the behavior of 
9s, gets completely modified by the presence of fermionic 
vacuum correction in the PQMVT and QMVT model 
and instead of dropping down to negative values (as 
in the PQM and QM model) in the chiral symmetry 
restored phase, 9s approaches the ideal mixing angle 
value smoothly similar to the temperature variation of 
the pseduoscalar 9p mixing angle for computations with 
the axial anomaly. 

In the chirally symmetric phase of the PQM and QM 
model in Fig 6(b) the scalar mixing angle drops down to 
9s ^ —51° (—54°) for higher temperatures in the presence 
(absence) of C/yi(l) axial symmetry breaking term. This 
drop happens around T/T^ - 1.5(1.9) in the PQM(QM) 
model for non zero c and similar drop for the calculations 
without axial anomaly occurs at a little higher value of 
T/T^. It is already reported and discussed in Ref . [28l. [soj 
that this behavior results because the masses of the 
physical a and /o anti cross and the nonstrange - strange 
(fivs — as) system masses cross in the close vicinity 
of the above mentioned reduced temperatures. After 
anti crossing the physical a becomes identical with pure 
strange quark system as while the physical /o becomes 
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FIG. 7: Shows the mass variations for the physical r;', rj and the nonstrange-strange rjMS, complex, on the reduced 
temperature scale (T/T^) at zero chemical potential (/i = 0). Fig |8(a)| shows the results for PQMVT and QMVT model 



and line types for mass variations are labeled, 
line types. 



Fig 8(b) shows the mass variations for the PQM and QM model with labeled 




FIG. 8: Shows the mass variations for the physical a, fo and the nonstrange-strange crjvs, as complex, on the reduced 
temperature scale {T/T}) at zero chemical potential (/i = 0). Fig 8(a) shows the results for PQMVT and QMVT model and 
line types for mass variations are labeled. Fig 8(b) shows the mass variations for the PQM and QM model with labeled line 
types. Here, the masses of the physical a and fo anti cross and the nonstrange-strange apfs — as system masses cross. 



degenerate with the pure nonstrange quark system (tns- 
In order to show this crossing-anticrossing behavior 
in the presence of axial anomaly, we have plotted in 
Fig |8(b)[ the PQM and QM model mass variations 
for the physical a and /o and the nonstrange-strange 
CATS) complex. Since the effect of fermionic vacuum 
fluctuation drastically modifies the 9s behavior for higher 
temperatures, the masses of the physical a and /o do not 
anti cross and the nonstrange - strange {uns~'^s) system 
masses do not cross for higher values on the reduced 
temperature scale and the a becomes identical with the 
pure nonstrange quark system uns while the physical /o 
becomes degenerate with the pure strange quark system 



CTs, in the PQMVT and QMVT model plots in Fig 8(a) 



VI. SUMMARY AND CONCLUSION 



In the present work, we have investigated how the 
inclusion of properly renormalized fermionic vacuum 
fluctuation in the 2-1-1 flavor QM and PQM models, 
modifies the finite temperature behavior of masses and 
mixing angles of scalar and pseudoscalar mesons. It 
has been explicitly shown that expressions for the 
model parameters, meson masses and mixing angles, 
do not depend on any arbitrary renormalization scale. 
We explored the qualitative and quantitative effects 
of fermionic vacuum correction, on the emerging mass 
degeneration patterns in the temperature variations of 
masses of the chiral partners in pseudoscalar (tt, 77, rj' , 
K) and scalar (cr, oq, fo,n) meson nonets. From the mass 
convergence patterns, we identified chiral symmetry and 
C/^(l) restoration trends and compared them in different 
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model scenarios. 

The temperature variation of the nonstrange, strange 
condensates and Polyakov loop field $ at /i = has 
been obtained from the gap equation in all the models 
QM, QMVT and PQM, PQMVT. Comparisons indicate 
that the temperature of the chiral crossover transition 
in the nonstrange sector shifts to the higher values in 
the influence of fermionic vacuum correction. The 
in the QMVT model increases by 25 MeV over its QM 
model value. We know that the Polyakov loop potential 
delays the chiral transition and increases significantly 
in the PQM model. Further one gets coincidence of 
the deconfincment crossover transition with the chiral 
crossover transition in the non-strange sector in the PQM 
model. We find that though the effect of fermionic 
vacuum fluctuation does not lead to any change in the 
deconfincment crossover transition temperature (T* = 
205.6 MeV in both the models; PQM and PQMVT), its 
influence spoils the above mentioned coincidence of the 
chiral transition with the deconfincment transition and 
registering an increase of 11 MeV over the T*=205.6 
MeV, the in the PQMVT model becomes 216.5 
MeV. The ax variation which is sharpest in the T/T'^ 
= 0.8 to 1.4 interval in the PQM model becomes 
smoother in the PQMVT model and in the absence of 
Polyakov loop, the QM model ax temperature variation 
becomes a lot more smooth in the QMVT model only 
because of the fermionic vacuum correction. It is 
already known that the enrichment of QM model with 
Polyakov loop potential, leads to a larger and significant 
melting of the strange condensate [23|- The melting 
of the ay is further enhanced due to the presence of 
fermionic vacuum correction and we obtain the largest 
but smoother melting of strange condensate in the 
PQMVT model. The interesting physical consequence of 
the earlier and larger but smoother melting of the strange 
condensate will be an early setting up of a smoother mass 
degeneration trend in the masses of the chiral partners 
[K^ k) and [rj, fo) and an early emergence of a smoother 
Ua{^) restoration trend. 

The sharpest mass degeneration trend obtained in the 
PQM model computations for (cr,7r) and (ap,?]') mesons, 
becomes quite smooth in the PQMVT model and the 
most smooth mass degeneration is seen in the QMVT 
model due to the absence of Polyakov loop potential. 
We thus conclude from the behavior of these chiral 
partners that the net effect of the inclusion of fermionic 
vacuum correction in the PQVI and QM model, is to 
make a smoother but more effective occurrence of chiral 
5t/L(2) X SUii(2) symmetry restoration transition in 
the nonstrange sector. Further in comparison to PQM 
and QM models, the degenerated masses of a and tt 
mesons in the PQMVT and QMVT models, show closer, 
smoother and more effective convergence towards the 
degenerate masses of oq and r]' for higher T/T^ > 1 
This behavior is a consequence of the largest but 
smoother and earlier melting of the strange condensate in 
the PQMVT model because the Ua{^) breaking anomaly 



effect that leads to the mass gap between the two sets of 
the chiral partners, (ct, tt) and (oq, 77') i.e. = ma- < 
TUao ~ ''^v' ^'^^ T /T^ > 1, is proportional to the strange 
condensate ay. Thus the incorporation of fermionic 
vacuum correction in the PQM and QM models also 
effects an early set up of Ua{^) restoration trend on the 
reduced temperature scale. 

Similar trend of smoother mass degeneration is 
noticed in the PQMVT and QMVT model temperature 
variations of the masses of the chiral partners (?7,/o) and 
(i^, k). Wc point out that in the QM and PQM models, 
the /o mass intersects the mass nearly at T/T^ = 
1.4(1.8) and becomes smaller than the rrir^ developing 
a kink like structure after crossing it for higher values 
of the reduced temperature T/T^ > 1.4(1.8). Further 
the fo meson becomes degenerate with [K, k) and 77 
again at T/T^ > 1.8(2.3). It is important to emphasize 
that this kink-like mass crossing behavior of /o meson 
altogether disappears from the PQMVT (QMVT) model 
mass variations of fo and 77 mesons due to the noteworthy 
effect of the fermionic vacuum correction. Here the fo 
meson becomes degenerate with the (K, k) and rj mesons 
earlier at T/T^ > 1.6(1.9) and remains so forever. This 
trend of mass degeneration reflects the effect of fermionic 
vacuum correction on the chiral symmetry restoration in 
the strange sector and it results again due to the smooth 
(smoother) but largest (larger) melting of the strange 
condensate in the PQMVT(QMVT) model. 

For T/T^ > 1 in the presence of axial anomaly, the 
pseudoscalar mixing angle Op approaches its ideal value 
more smoothly in the PQMVT and QMVT model when 
compared with the corresponding result in the PQM and 
QM model. The smallest mass difference between the 
771^' (which is degenerate with Tn^^g) and rri,, (which is 
degenerate with m^.,) for T/T^x > 1 results in PQMVT 
model. This again is the evidence of the setting up of the 
most effective ?7a(1) restoration trend in the PQMVT 
model for the pseudoscalar sector. 

Due to the effect of fermionic vacuum correction, the 
scalar mixing angle 0s value in vacuum (i.e. at T = 0) 
decreases to 11.98(14.75) degree in the presence (absence) 
of axial anomaly in the PQMVT and QMVT models from 
its value of 19.86(21.5) degree for PQM and QM models. 

In the chirally symmetric phase of the PQM and QM 
model, the scalar mixing angle drops down to 9s ^ 
—51° (—54°) for higher temperatures in the presence 
(absence) of J7a(1) axial anomaly. Due to this drop, 
the masses of the physical a and fo anti cross and the 
nonstrange - strange {a^s ^ <^s) system masses cross in 
the QM and PQM model. The most striking effect of 
the fermionic vacuum correction is seen in the complete 
modification of the temperature variation of the scalar 
mixing angle 9s in the chirally symmetric phase of the 
PQMVT and QMVT model. Instead of dropping down 
to large negative values, it smoothly approaches the ideal 
value of {9s ^ 35°) mixing angel for higher temperatures 
in the presence as well as absence of axial anomaly. As 
a consequence, masses of the physical a and fo do not 



16 



anti cross and the nonstrangc-strange (cfns ~ f s) system 
masses do not cross for higher values on the reduced 
temperature scale in the PQMVT(QMVT) model, a 
becomes identical with pure nonstrange quark system 
aNS while the physical /o becomes degenerate with the 
pure strange quark system as- 
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